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Abstract 

In connection with the semi-inclusive polarized DIS, it is proposed to consider 
the first Mellin moments Aq of the polarized quark and antiquark densities, instead 
of the respective variables 6q(x), local in Bjorken x themselves. This gives rise to a 
very essential simplification of the next to leading order (NLO) QCD and, besides, 
allows one to use the respective QCD sum rules. An expression for Au — Ad in NLO 
is obtained which is just a simple combination of the directly measured asymmetries 
and of the quantities taken from the unpolarized data. 



1 Delivered on October 10 at Physics Workshop "Compass week in Dubna" (JINR, Dubna, 2000) 
under the title "Polarized sea-quark flavor asymmetries and COMPASS." 
2 E-mail address: shevch@nusun.jinr.ru 
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Investigation of the quark structure of the nucleon is one of most important tasks of 
modern high energy physics. In this respect deep inelastic scattering (DIS) is of special 
importance. Thus, the very impressive result of the New Muon Collaboration (NMC) 
experiment was obtained in 1991, when the unpolarized structure functions of the proton 
and neutron, F% (x) and F£(x), were precisely measured within a wide range of Bjorken's 
x, and, it was established that the integral f c \ —[F%{x) — F£(x)\ does not equal 1/3 
(Gottfried sum rule) but has a much smaller value 0.235 ± 0.0026. This means that the 
densities of u and d sea quarks, u(x) and d(x), in the proton have different values, and 

f 1 dx \d(x) - u(x)] = 0.147 ± 0.039 ^ 0. 
Jo 

In polarized DIS, instead of the unpolarized total q = q^ + q±, sea q and valence 
qy = q — q quark densities, the set of the respective polarized quantities Sq(x, Q 2 ) = 
q^(x,Q 2 ) - qi (x,Q 2 ), 5q(x,Q 2 ) = q\x,Q 2 ) - q^x.Q 2 ) and 5q v (x,Q 2 ) = 5q(x,Q 2 ) - 
Sq(x, Q 2 ) is the subject of the investigation. So, the question arises: does the difference 
between the polarized u and d sea quark densities 8u(x, Q 2 ) — Sd(x, Q 2 ) also differ from 
zero? Recently, a series of theoretical papers appeared ([1-4]) where it was predicted that 
the quantity 5u(x, Q 2 ) — 8d(x, Q 2 ) does not equal zero. However, the model-dependent 
results for Su(x, Q 2 ) — 8d(x, Q 2 ) essentially differ each from other in these papers. So, it 
is very desirable to find a reliable way to extract this quantity directly from experiment 
data. For this purpose it is not sufficient to use just the inclusive polarized DIS data, and 
one has to investigate semi-inclusive polarized DIS processes like 

p. + p(d) -> p + h + X. 

Such processes provide direct access to the individual polarized quark and antiquark 
distributions via measurements of the respective spin asymmetries. f\ 

Unfortunately, the description of semi-inclusive DIS processes turns out to be much 
more complicated in comparison with the traditional inclusive polarized DIS. First, the 
fragmentation functions are involved, for which no quite reliable information is available^. 
Second (and this is the most serious problem), the consideration even of the next to 
leading (NLO) QCD order turns out to be extremely difficult, since it involves double 
convolution products. So, to achieve a reliable description it is very desirable, on the one 
hand, to exclude from consideration the fragmentation functions, whenever possible, and, 
on the other hand (and this is the main task), to try to simplify the NLO consideration 
as much as possible, without which one can say nothing about the reliability and stability 
of results obtained within the quark-parton model (QPM). 

It is well known that within QPM one can completely exclude the fragmentation 
functions from the expressions for the valence quark polarized distributions 5qy through 
experimentally measured asymmetries. To this end, instead of the usual virtual photon 
asymmetry A^ N = A\ N (which is expressed in terms of the directly measured asymmetry 
Ae X p = (n^ — n^)/(n^ + n^) as A\ N = (P B PTfF))~ 1 A^ xp ), one has to measure so called 

"difference asymmetry" A^~~ h [6] (see also [5,7]) which is expressed in terms of the 

3 Such a kind of measurements were performed by SMC and HERMES experiments and are also 
planned by the COMPASS collaboration. 

4 For discussion of this subject see, for example [5] and references therein. 
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respective counting rates as 



Ah-h ( n 2. \ = 1 (n^ n^) (wfo. n^ f 
N 1 ' W ' J P^Pt/^D - n\) + (n* - 4 T ) ' 1 J 

where the event densities ^jm) = dN^^^/dz, i.e. n^i"^^ are the numbers of events 
for antiparallel (parallel) orientations of here muon and target nuclear (proton or deutron 
here) spins for the hadrons of type h registered in the interval dz. Coefficients Pb and 
Pt, f and D are the beam and target polarizations, dilution and depolarization factors, 
respectively, (for details on these coefficients see, for example, [8,9] and references therein). 
Then, the QPM expressions for the difference asymmetries look like (see, for example, 
COMPASS project [10], appendix A) 
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A p 


— 7T 


4Suy — Sdy 
Au v — d v 


ATT+ 


— 7T~ 


buy + Sdy 
u v + d v 


aK+ 


-K- 


5u v _ 

Uy ' 


A K+ - 

A d 


-K- 


A 7T"*" — 7T — 

— A d J 



45dy — Suy 
Ady — Uy 



(2) 

i.e., on the one hand, they contain only valence quark polarized densities, and, on the 
other hand, have the remarkable property to be free of any fragmentation functions. 

All this is very good, but we are interested here in the sea quark polarized distributions, 
and, besides, the main question arises - what will happen with all this beauty in the next 
to leading order QCD? 

We propose to investigate the integral quantities, namely, the first Mellin moments 
M l (8q) = Jq 1 dx [5q(x)] = Aq (q = u,d, s, ...) instead of the local polarized quark densities 
8q(x) themselves. This provides very essential advantages: 
First . 

Even if the local quantity has a very small f] value at each point x, the integral of this 
quantity over the whole range of x-variables may already have quite a considerable value, 
and, one can hope that QPM turns out to be a good approximation for integral quantities 
like 

Au-Ad=f dx \Su(x) - 5d(x)\. (3) 
Jo 

An argument in favor of such a hope (for (3)) is the circumstance that all the model 
predictions [1-4] have one common feature: the local quantity 5u(x) — 5d(x) does not 
change sign when x varies over its entire range < x < 1. 
Second . 

To investigate integral quantities like (3) one can use QCD sum rules. In particular, one 
can apply such a well established sum rule as the Bjorken sum rulef\ 

1 dx[gf - 9l ] = l 9 -A {1 _ + Q{ 2))j (4) 

o 6 gy 7i 



5 Notice, however, that the latest theoretical paper [4] on this subject predicts that the difference 
between the polarized densities Su and Sd should be even more significant than the difference between 
the unpolarized sea quark densities: \5u — 8d\ > \u — d\. 

6 Throughout the paper, all the quantities considered in NLO are given in the MS scheme. 
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g A / gv = 1.2537 ±0.0028 

to express the quantity Au — Ad of interest via the quantity Auy — Ad v which, in turn, 
is expressed via the measured difference asymmetries A^ + ~ n and _7r . 

Third (and we consider this the most important advantage of the proposed procedure) 
Application of the Mellin moments, instead of the local quantities themselves, results in a 
remarkable simplification of the NLO QCD consideration of the semi-inclusive polarized 
DIS, that is extremely complicated in terms of the local quantities. 

Thus, let us consider the NLO [11] expression for the structure function g\ 

Q 2 ) = I E < (sq + ^pfo Sq + c ° Sg] ) {x > g2) ' (5) 



where 



(c®/) W ./;^ c (£)/(,) ( 6) 

is the definition of the convolution product. From now on we will use the well known 
remarkable property of the Mellin n-th moments 

M n (f)= ( l dxx n - l f(x) (7) 
Jo 

to split the convolution product (6) into a simple product of the Mellin moments of the 
respective functions: 



M n [C ® /] = f 1 dxx n ~ l f (-) f(y) = M r 
Jo Jx y \yj 



{C)M n {f). (8) 



So, taking the first Mellin moment of Eq. (5) and using the expressions for the Mellin 
moments of the respective NLO (MS) Wilson coefficients 

M 1 (C 9 ) = -2, M\C g )=0 

one obtains [11] in NLO QCD: 

M 1 \g[] ee jf 1 dx g>(x, Q 2 ) = l -Y, ^ (l " [ dx Sq (9) 

and the same for g™ with the substitution u <-> d. Substituting the last expressions for 
M 1 ^] and M 1 ^™] into the Bjorken sum rule (4), one can see that the a s dependent 
multipliers (1 — a s (Q 2 )/n) cancel out precisely in the left- and right-hand sides, and, one 
arrives at the simple relation between the polarized densities of sea and valence quarks 

/ dx (5u — 5d) — — — [ dx (5uv — Sdy), (10) 

Jo 2 gy 2 jo 

or, in the notation used here, 

Au- Ad= - —-\ (Au v - Ad v ). (11) 
2 gv 2 
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Thus, the relation between the first Mellin moments of the polarized sea and valence 
quark distributions has a very simple form and does not contain a s dependence at all (i.e. 
is an exact relation at least up to 0(a 2 ) corrections). 

With such a simple relation between Au — Ad and Auy — Ady at our disposal, the 
next step is to establish the relation between the Mellin moments Auy and Ady and the 
experimentally observable difference asymmetries in NLO QCD. For this purpose, 

one can use the following relations [5,12,13] for the difference asymmetries 

N/h _ N/h 

Afr\x,<?;z)= 9 * (N = p,n,d), (12) 

t 1 - b x 

where the semi-inclusive analogs of the structure functions and Ff, functions g^ h and 
F^ h , are related to the respective polarized and unpolarized semi-inclusive differential 
cross-sections as follows [12] 

dxdydz dxdydz Q 2 ^ v) 9i ( i iQ )i ( ) 

dV N Wi + o-jo^^ ^ 



dxdydz Q 2 



where 



pN/h _ F N/h 1 -V F N/h ,. -, , 

*1 = *1 + l + (l_j,)2 ■ ^ 



The semi-inclusive structure functions g^ n ^ h are given in NLO by 



q,q Z7F 

+ (L^)®^c gq ®D h g 

q,9 Z7F 

+ 5g®^C qg ®C£e 2 q D h q l (16) 



9,9 



n/h p/h 

9i = 9i 



(17) 



where 



[A*B<aC\{x,z) S / I^A(^)B( X ',Z')C(^) (18) 

V 

is the double convolution product. The respective expressions for 2F^ n ^ h have the same 
form with the substitution 5q — > q, 5C — > C. The expressions for the Wilson coefficients 

SC qq(qg,gq) aIld Cqq{qg,gq) = C \ q {qg, gq ) + 2 [( X " V) I ' i 1 + i 1 ~ Vf)\ C ,gq) CaI1 be foUIld ' for 
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example, in [12], Appendix C. 

It is remarkable that due to the properties of the fragmentation functions: 

D x ee Df = Df =Df =Df, 

D 2 = Df = Df =DZ~ =Dl\ (19) 

in the differences g\°^ — and Ff^ — fif (and, therefore, in the asymmetries 

Ap + ~ w and A 7 ^^ ) only the contributions containing the Wilson coefficients 5C qq and 
C qq survive. However, even then the system of double integral equations 

A n+-n- (rr n 2. x (4<W - 5dy)[l + ®a 8 /(2'K)6C qq ®](Di - D 2 ) 2 . 

p [x ^> z) - (4u v - d v )[l + ®a s /{2-K)C qq ®]{D x - D 2 ) ^ ^ 
Al + --(x,Q 2 ;z) = A; + --(x,Q 2 ;z)\ Uv ^ dv (20) 

proposed by E. Christova and E. Leader [5], is extremely difficult to solve with respect to 
the local quantities 5uv(x,Q 2 ) and 5dy(x, Q 2 ). Besides, the range of integration V used 
in ref. [5] has a very complicated form, namely: 

x 



— < x' < 1 with z < z' < 1, (21) 
x + (1 — x)z 

if x + (1 — x)z >1, and, additionally, range 

x < x' < xj[x + (1 — x))z 

with x(l — x')/ (x'(l — x)) < z' < 1 if x + (1 — x),2 < 1. So, one can see that here even 
application of the Mellin moments cannot simplify the situation. 

Such enormous complification of the convolution integral range occurs if one in- 
troduces (to take into account the target fragmentation contributions 7 and to exclude 
the cross-section singularity problem at Zh = 0) a new hadron kinematical variable 
z = Eh/E^(l — x) (7P cm. frame) instead of the usual semi-inclusive variable Zh = 
(Ph)/(Pq) = (lab. system) Ej l /E 1 . However, both problems compelling us to introduce 
z, instead of Zh, can be avoided (see, for example [12,13]) if one, just to neglect the target 
fragmentation, applies a proper kinematical cut Z < Zh < 1, i.e. properly restricts the 
kinematical region covered by the final state hadrons^. Then, one can safely use, instead 
of z, the usual variable z^, which at once makes the integration range T> in the double 
convolution product (18) very simple:0 x < x' < 1, z^ < z' < 1. Note that in applying 
the kinematical cut it is much more convenient to deal with the total numbers of events 

Nf l(n) (x,Q 2 ) = f dz h n h n{m {x,Q 2 -z h ) (22) 



within the entire interval z < z-^ < \ and the respective integral difference asymmetries 

1 (JVft - Nfr) - (Mft - Af T ) 



A h N h (x, Q A 



PbPtJD (N^ - JV* ) + (JVft - 



z 



7 Then, one should also add the target fragmentation contributions to the right-hand side of (16). 
8 This is just what was done in the HERMES and COMPASS experiments, where the applied kine- 
matical cut was Zh > Z = 0.2. 

9 Namely the such range was used in the early seminal papers [14] (see also [12]). 

10 Namely the integral spin symmetries A^ N = f^, dz} x g^ h j f% dz^ p^ h were measured by SMC and 

HERMES experiments (see [8,9] and also [13]). 
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Iz dz h(9i 



N/h 



N/h, 

9i ) 



Jhdz h (F»"*-Fl 



\N/h, 



(N = p,n,d), 



(23) 



than with the local in Zh quantities n^u^)(x, Q ; Zh) and A N (x, Q ; z/i). So, the expres- 
sions for the proton and deutron integral difference asymmetries assume the formF] 



(x,Q 2 



{A8u v - Sd v ) Si dz h [l + ^SC^jD, - D 2 
{4u v - d v ) fedz h [\ + ®%;C qq ®}(Di - D 2 ) 



(24) 



A 77 



(x,Q 2 



-1 j.. \f jp/k + 



Iz dz h Mi 



p/tt 

9i 



(9i 



n/ir + 



n/ir 

9i 



)] 



Hdz h [{F^-F{ 



n/ir~^ 



F^)} 



_ (Su v + 5d v ) j l z dz h [l + ®7^5C qq ®}{Di - D 2 

{Uy + dy) Jldz h [l + ®^C qq ®]{D l - A>) 

and the double convolution reads 



[A ® B <g> C] 



1 dx' r 1 dz' ^fx 

X 1 Jz h Z 1 \x' 



B(x',z')C( Z -±). 



(25) 



(26) 



Now, application of the first Mellin moment to the difference asymmetries n (x, Q 2 )\z 
and A^ + ~ n (x, Q 2 )\z, given by (24) - (26), becomes extremely useful and allows one to 
obtain a system of two purely algebraic equations for Auy = Jq 1 dx 5uy and Ad v = 
Jq dx 5dy: 



(4Auy - Ad v ) (Mi - M 2 



a exp 



(27) 



(A U y + Ady)(M 1 - M 2 ) = A, 



exp 



:28) 



with the solution 



lA^ + A e d xp v , 
Au v = --^ Ad x 



5 Mi - M 2 

Here we introduce the notation 

-l 



1 AAT - A e p x P 
5 M 1 - M 2 



a exp 



jexp 



dx Ap 



and 



Mi 



(Auy-dy) dz h [l + ®-±C qq ®](D 1 -D 2 ), 

JZ 2,71 
rl OL ~ 

{uy + dy) J z dz h [l + ® 1 ^C qq ®}{D l - D 2 ), 



Ml 



Mo = M 



d 



Ml 



Ml , 
MJ + , 



(29) 



(30) 



(31) 



11 To obtain (25) one uses that the deutron cross section is the sum of the proton and neutron cross 
sections, which is valid up to corrections of order O(wd), where ujd = 0.05 ± 0.01 is the probability to 
find deutron in the _D-statc. 
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where 



M h q (Q 2 ) = f dz h \D h q (z h} g 2 ) + ^f^r &C qq (z')D h q (^, Q 

J Z ^7T J Zh Z Z 



(32) 



with the coefficient 



AC qq (z) = / dxSCgg(x,z), (33) 



o 



that is given in Appendix. Thus, using the relation (11) between Am — Ac? and Auy — Ady 
one gets, eventually, a simple expression for Am — Ad = Jq dx (5u(x, Q 2 ) — 5d(x, Q 2 )) in 
terms of experimentally measured quantities, that is valid in NLO QCD : 

Au — Ad = - — — -—j (34) 

2 g v 10(M 1 -M 2 ) v ' 

It is easy to see that all the quantities present in the right-hand side, with the exception 
of the two difference asymmetries A 7 ^' 1 * \z and A\ ~ n \z (entering into Ap Xp and Atp ', 
respectively) can be extracted from unpolarized\^\ semi - inclusive data and can, thus, 
be considered here as a known input. So, the only quantities that have to be measured 
in polarized semi- inclusive DIS are the difference asymmetries A^ + ~ n \z and A\ _7r \z 
which, in turn, are just simple combinations of the directly measured counting rates. 

In conclusion, we would like to stress that application of the Mellin moments, instead 
of the local polarized densities, happens to be very fruitful not only in the case of light u- 
and d-quarks, but also for investigation of polarized strangeness in the nucleon (a paper 
is now in preparation). Besides, we also plan to apply this procedure to the transverse 
asymmetries in the nearest future. 

At present, a proposal for measurement of Au — Ad , based on the above described 
procedure, is being prepared for the experiment COMPASS in collaboration with the 
group of INFN - sezione di Torino and of Dipartimento di fisica generale "A.Avogadro" 
of the Torino University. 
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Appendix: Mellin moments of polarized semi-inclusive DIS Wilson 

coefficients. 

The NLO (MS) coefficient SC qq has the form (see [12], Appendix C) 

5C qq = C 1 qq -2C F (l-x)(l-z) (C F = 4/3), (35) 

where 

12 With the standard and well established assumption that the fragmentation functions do not depend 
on the spin. Then, the unpolarized fragmentation functions D can be taken either from independent 
measurements of e + e~ - annihilation into hadrons [15] or in hadron production in unpolarized DIS [16] 
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CL = C F \-SS(l-x)S(l-z)+S(l-x)\P qq (z)\n^ T + L 1 (z) + L 2 (z) 



M F 



+ (1 -z)\ + 5(1 -z) 
1 1 



Q 2 



+ 



+ 2 



1 + z 



1 + x 



(l- x ) + (l-z) + (l-x) + (1-z). 



+ 2{l + xz) \ , 



(36) 



L ± (z) = (1 + z 2 ) 



'ln(l -z)\ 



1 + z2 

(T^jT + 2 5(1 -^ 



, L 2 {z) = 



1 + z 2 



1 — z J + 1 — z 

and the " +" prescription is given by 

fdz f(z)(g(z)) + = fdz [f(z) - f(l)} g(z) 
Jo Jo 

Simple calculation of the n-th moment M n (SC qq ) gives 



lnz, 



(37) 



(38) 



M n (8C qq ) = C F {5(l-z) 



+7^ 6^ + 3 



3, Q 2 , Q 2 /3 n 1 1 
-8 + - In + In - - 27 

2 Mp M 2 V2 f n 1 + n 



- 2* (n) 
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v n 2 (l + n) 2 7 
+ 3^ 2 (n + 2) -6 



+ 67 ( - + + 7T 2 + 12 7 * (n) + 3^ 2 (n) 

n 1 + nj 

d^(n)\ \ a/o o \ 1 



n(n + 1) 



-[7 + *(n)] + (1 + z)[7 + *(n)] 



(1-^)+ 

1 /l 1 \ (1 : 
' +2(- + 



(1 — z) + \n 1 + n 



n 1 + n 



+ P„Wh.g + L.W + + (1 - -(^^^ 



where = r(z)/r(z); 7 ~ 0.577216 is the Euler constant. 
For example, the first two moments are 



M x (SC gq ) = AC 9? = C F 



l + 2z 



2(l-z) + 



re 2 3, g_ 2 







(39) 



M 2 (<SC OT ) = C F 

q 2 , n g 2 



' + 2z - — - - + h(l - z)( - 41 + 2vr 2 - 



6 (1-z) 

Q 2 



8ln M 2+ 9 W) + Pqq{z) ln W + Ll{z) + U{Z) + 3 {1 ~ Z) 



(40) 
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